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摘要 :目的 描绘 小 鼠 雌 性 生殖 干细胞 (FEGSCs) 和 胚胎 干细胞 (ESCs) 全 基因 组 范围 内 选择 性 多 聚 腺 昔 酸 化 (APA) 位 点 ,通过 细 
胞 间 的 比较 ,分 析 生 殖 干 细胞 特异 性 的 APA 位 点 对 其 生物 学 特性 的 影响 。 方 法 利用 本 实验 室 之 前 建立 的 基于 高 通 量 测序 的 转 
录 本 3 末端 poly(A) 位 点 捕获 方法 3Tseq 确 定 并 比较 两 个 细胞 系 中 所 有 APA 位 点 ,通过 DAVID 对 具有 APA 位 点 差异 的 基因 进 
行 Gene Ontology 分析。 结果 在 两 种 细胞 系 中 共 确 定 16 973 个 基因 的 50 243 个 APA 位 点 ,发 现 FGSCs 中 相 较 于 ESCs 3'UTR 
发 生长 度 变 化 的 基因 1148 P , Horn 795 47 (6696)3' UTR 26/8 ,353 4- (3496 ) ELS] 3' UTR 变 长 ,与 生殖 发 育 密切 相关 的 基因 3” 
UTR 存 在 显著 的 缩短 现象 。 结 论 FGSCs 存 在 与 ESCs 不 同 的 APA 谱 ,APA 介 导 的 3”UTR 变化 参与 生殖 发 育 相关 的 生物 学 过 


程 ,揭示 了 生殖 干细胞 的 基因 转录 后 调控 机 制 。 
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Abstract: Objective To perform a genome-wide alternative polyadenylation (APA) profiling in both mouse female germline 


stem cells (FGSCs) and embryonic stem cells (ESCs) and explore the role of germline-specific APA in the biological behaviors 
of FGSCs. Methods We used a high-throughput sequencing-based method 3T-Seq to profile the genome-wide 3' termini of the 
transcripts and delineate all the APA sites in mouse FGSCs and ESCs. The genes with altered APA sites in FGSCs compared 
with ESCs were analyzed with DAVID Gene Ontology tool for their biological roles. Results We identified a total of 50243 APA 
sites in 16973 genes. In FGSCs, 1148 genes were shown to have alterations in 3'UTR length, among which 795 ( 66%) genes had 
shortened and 353 (34% ) had lengthened 3'UTR. Some of the genes with shortened 3UTR were involved in germ cell 
development. Conclusions Our genome-wide APA profiling analysis reveals a cell type-specific APA alternation in FGSCs, and 


APA-mediated 3'UTR alteration contributes to germline-related biological process. This study provides a framework for 
understanding the post-transcriptional regulation mechanisms in FGSCs. 
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选择 性 多 聚 腺 背 化 (APA) 是 一 类 很 重要 的 基因 转 
录 后 表达 调控 方式 ,通过 调节 转录 本 3 " 非 翻译 区 (3” 
untranslated region, UTR) 的 长 度 来 影响 包括 mRNA 的 
稳定 性 、 亚 细胞 定位 和 翻译 效率 等 在 内 的 多 个 生物 学 过 
TE, APA 在 生物 体 生 长 发 育 和 疾病 发 生发 展 过 程 中 
受到 高 度 有 序 的 调控 。 在 发 育 生 物 学 领域 ,全 基因 组 
范围 内 的 APA 谱 表现 出 细胞 类 型 特异 性 :在 多 能 性 的 
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豚 胎 干细胞 中 ,基因 倾向 于 选取 近 端 的 多 聚 腺 苷 酸 位 点 
(PAS) 而 产生 更 短 的 3"UTR 来 维持 干细胞 特性 ;而 终端 
分 化 细胞 更 多 取 用 远 端 PAS "7; 由 APA 调控 的 3’'UTR 
的 长 短 也 呈现 发 育 阶段 性 :小 鼠 胚 胎 发 育 伴随 着 全 局 3" 
UTR 的 逐步 增长 ; 果 晶 胚胎 发 育 过 程 中 的 3”"UTR 变 化 
趋势 与 之 相似 外 ;而 斑马 鱼 的 胚胎 发 育 中 mRNA 的 3 
UTR 呈现 “V"” 字 形 的 动态 变化 ,在 合子 形成 后 基因 组 被 
活化 表达 的 短 时 间 内 3”UTR 缩短 ,之 后 迅速 变 长 ,日 伴 
随 着 3”UTR 多 样 性 的 由 丰 到 简 ”。 

小 鼠 肉 性 生殖 干细胞 (mFGSCs) 是 一 类 单 能 性 成 体 
二 细胞 ,能 在 保持 自我 更 新 则 时 定向 分 化 为 卵细胞 “? ， 
是 生殖 医学 和 干细胞 研究 领域 的 重要 研究 对 象 。 虽 然 
从 小 鼠 体 内 成 功 分离 并 建立 细胞 系 ”, 但 FGSCs 的 分 


201712.02111v1 


chinaXiv 


* 158 - J South Med Univ, 2016, 36(2): 157-162 


ChinaXiv& f'ERBTII 


http://www.j-smu.com 


子 生物 学 研究 目前 仍 处 于 起 步 阶段 ,其 分 子 水 平 上 的 调 
控 机 理 至 今 尚 不 清晰 。APA 所 产生 的 不 同 3”"UTR 作 为 
一 种 重要 的 表 观 遗传 调控 手段 在 生殖 干细胞 中 所 发 挥 
的 作用 也 未 被 深入 研究 。 本 文 应 用 本 实验 室 建立 的 基 
于 高 通 量 测序 的 3? 端 poly(A) 位 点 捕获 技术 3T-seq'a ， 
在 全 基因 组 范围 内 分 别 对 小 鼠 肉 性 生殖 干细胞 和 豚 胎 
干细胞 (ESCs)APA 谱 图 进行 了 分 析 比 较 , 以 探讨 该 转 
录 后 修饰 对 FGSCs 特 性 的 影响 ,为 进一步 丰富 和 完善 
生殖 干细胞 分 子 生 物 学 研究 提供 新 的 内 容 。 


1 材料 和 方法 
1.1 细胞 培养 

小 鼠 肉 性 生殖 干细胞 (由 上 海 交 通 大 学 吴 际 教授 惠 
赠 ), 小 鼠 胚 胎 干 细胞 敌 均 培养 于 用 0.1% 明 胶 包 被 的 6 筷 
板 培养 四 (Corning),FGSC 培 养 基 为 MEM-a(Gibco)， 
含有 终 浓 度 为 10% 胎 牛 血清 (Gibco) ,1 mmol/L 丙酮 酸 
钠 ,1 mmol/L 非 必需 氨基 酸 (Gibco),1 mmol/L L-Zi24 
酰胺 ,0.1 mmol/L B- 琉 基 乙 醇 ,20 ng/mL LIF(Peprotech) , 
10 ng/mL 的 EGF (Peprotech) , 10. ng/mL GDNF 
(Peprotech) , 10 ng/mL bFGF (Peprotech) ; ESC 培养 基 
为 :高 糖 型 DMEM(Gibco) 15% 胎 牛 血清 .100 U/mL $$ 
HEREZ (Gibco) .1000 U/mL LIF(Chemicon) ,0.1 mmol/L 
B- 琉 基 乙 醇 1% 非 必需 氨基 酸 ;新 鲜 配制 的 培养 基 用 
0.22 hm 的 滤器 (Millipore) 过 滤 后 4 CC 保存 。 细 胞 培养 
置 于 5% CO,,37 %C 培 养 箱 培养 ,每 2 d 更 换 培养 基 。 细 
胞 生长 至 80%~90% 密 度 时 ,以 1:2~3 的 比例 传代 ,经 3~ 
4 代 收 集 细胞 用 于 提取 总 RNA。 
1.2 细胞 总 RNA 的 提取 

将 生长 状况 良好 的 细胞 弃 去 培养 基 , 用 PBS 洗 1 
次 ,每 个 6 孔 板 加 入 1 mL Trizol regent(Invitrogen) 收 集 
细胞 ,反复 吹 打 使 贴 壁 细胞 破碎 溶 于 Trizol 后 转移 至 新 
的 RNase-free 的 1.5 mL 离心 管 中 , 用 于 立即 抽 提 总 
RNA SÉ T -80 保存 。 总 RNA 抽 提 参 照 试 剂 盒 说 
HE, RNA YLE F RNase Free 7K (Invitrogen) ,用 
Nanodrop 2000(Thermo) 定 量 并 在 1.5% 琼 脂 糖 凝 胶 中 
检验 RNA 完 整 度 。 
1.3 3T-seq 文 库 构 建 

文库 构建 参照 本 实验 室 之 前 报道 的 方法 ” ,简单 来 
说 :50 hg 总 RNA 与 预先 包 被 了 带 生 物 素 (Biotin) 修 饰 
的 反 转 录 引 物 Gsul-Oligo(dT)x 的 链 霉 亲 和 素 磁 珠 
M-280(Invitrogen) 孵 育 , 在 Super Script Ill (Invitrogen) 
反 转 录 酶 的 作用 下 ,用 含有 5-methylated-dCTP 的 dNTP 
Mix 合成 第 1 链 cDNA,5m-dCTP 封闭 新 合成 的 一 链 
cDNA 内 部 Gsu 工 的 识别 和 切割 位 点 ;经 过 二 链 合 成 后 
在 磁 珠 上 将 双 链 cDNA 用 片段 化 酶 (NEB ) 随 机 打 断 为 
200-400 bp, 只 留 下 与 磁 珠 相连 的 近 3’ 端 的 DNA 标签 


片段 ,用 GsuI(Fermentas) 来 酶 切 释放 。 纯 化 后 将 37UTR 
标签 文库 末端 补 平和 3 ? 端 加 A, 连 接 Ilumina 公 司 的 P3/ 
P7 测 序 接头 ,高 保 真 酶 PCR 放 大 ,筛选 200~500 bp 大 小 
DNA 作为 测序 文库 经 Illumina Hiseq 2500 进 行 高 通 量 
测序 ,测序 数据 已 上 传 至 EMBL-EBI ArrayExpres 数 据 
库 (accession number: E-MTAB-4038)。 
1.4 数据 分 析 

数据 处 理 方法 参照 之 前 的 报道 中 ,将 原始 数据 经 过 
C++ 脚本 的 过 滤 和 筛选 得 到 的 可 用 有 效 的 数据 ,应 用 
bowtie2 软件 比 对 到 UCSC 的 小 鼠 参考 基因 组 mm10， 
APA 位 点 由 迭代 聚 类 方法 鉴定 。 两 种 细胞 系 之 间 的 
3' UTR 的 变化 用 线性 趋势 分 析 方 法 进行 比较 。FDR 由 
R 软件 估计 ,基因 功能 分 析 Gene Ontology Hi DAVID 
(https://david.ncifcrf.gov) 完 成 。 


2 结果 
2.1 FGSCs 和 ESCs 细 胞 内 基于 3T-seq 的 全 基因 组 范围 
APA 位 点 的 定位 

高 质量 的 总 RNA 提取 之 后 按照 我 们 最 近 报 道 的 方 
法 5 构建 3Tseq 文 库 ( 图 1A) 并 进行 高 通 量 测 序 。 对 测序 
数据 进行 生物 信息 学 处 理 , 两 个 样品 共 得 到 约 9.2 M 
reads {EX clean data, 其 中 7.7 M( 约 84% )reads 能 够 比 
对 到 基因 组 上 ,包含 3.88 M 特 异性 reads( 占 50% )。 基 
Poly(A) 位 点 下 游 20 nt 序列 车 出 现 连续 8 个 A 或 含 
有 12 个 或 以 上 的 A, 则 作为 Internal Priming WEH VI 
减少 对 真实 APA 位 点 鉴定 的 影响 ,最 终 过 滤 得 到 3.1 M 
有 效 reads, 这 些 reads 有 接近 89.7%( 曲 线 下 方 的 面积 
集中 在 已 注释 的 转录 终止 位 点 (Transcription Terminal 
site, TTS) 附 近 , 如 图 2A 所 示 。APA 位 点 鉴定 方法 根据 
之 前 报道 中 ,最 后 共 确 定 50 243 个 APA 位 点 :FGSCs 细 
胞 系 25 431 个 ,ESCs 细胞 系 24 812 个 。 在 全 基因 组 范 
围 内 位 点 如 图 1B 所 示 , 具 体 统计 数据 见 表 1。 在 APA 
位 点 上 游 20~30 nt 范围 内 搜索 以 AAUAAA 为 序列 的 
PAS 六 聚 体 ,如 图 2B 所 示 , 除 了 CDS 区 的 APA 位 点 外 
均 能 找到 该 motif, 表 明定 义 APA 位 点 的 高 效 性 。APA 
位 点 在 基因 组 中 的 具体 分 布 如 图 2C,42.7% 位 于 已 知 3 
UTR 区 ,29.7% 落 在 转录 终止 位 点 (TTS) ,剩余 的 分 别 比 
对 到 TTS 下 游 1000 nt:4.3% ,编码 区 (CDS):11.2% ,内 
含 子 区 (Intron) : 2.9% , 非 编码 区 (Noncoding gene) : 
1.8% 和 基因 间隔 (Intergenic):7.5%。 
2.2 FGSCs 与 ESCSs 之 间 3?UTR 的 差异 分 析 

为 了 方便 进行 比较 ,我 们 采用 线性 趋势 分 析 (Test 
of linear trend) 方 法 “对 含有 多 个 APA 位 点 的 基因 均一 
化 处 理 后 引入 3”UTR 长 度 指数 (3’UTR length index, 
ULT) 来 定量 指 征 同 一 个 基因 3’UTR 的 长 度 变化 , 若 
ULI 为 正 值 , 则 表示 FEGSC 细 胞 系 的 3"UTR 相对 于 ESC 
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图 1 FGSCs 和 ESCs 中 APA 位 点 的 研究 

Fig.1 Identification of APA sites in FGSCs and ESCs. A: Gel electrophoresis of the total RNA and 3T-seq library 
for NGS. Lane 1: FGSCs; Lane 2: ESCs. M: 100 bp DNA Ladder; B: A genomic view of APA sites defined by 
3T-seq in IGV genome browser. Blue track: ESCs; Yellow track: FGSCs. 
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表 1 高 通 量 测序 数据 处 理 和 比 对 结果 汇总 
Tab.1 Summary statistics of sequencing data from Ilumina 
Hiseq2500 


FGSCs ESCs Combined 
Total reads after trimming 5332614 395632] 9288935 
Reads Mapped to genome 4189582 3579438  À 7769020 
Reads Uniquely mapped 2215930 1672353 3888283 
to genome 
Reads Passed internal 1746257 1346373 3092630 
priming filter 
Identified Poly(A) sites 25 431 24 812 50 243 
Genes hit by poly(A) sites 8464 8509 16 973 


的 变 长 ,反之 变 短 。 结 果 有 1 148 EDI 3 UTR 长 度 发 
生 明 显 变化 ,其 中 ULI 为 负 795 个 (66% ) ,ULI 为 正 353 
个 (34% )。 整 体 的 变化 趋势 如 图 3A 所 示 。 拥 有 不 同 
APA 位 点 的 基因 个 数 代表 着 细胞 APA 的 多 样 性 ,图 2D 
展示 了 这 一 分 布 情况 。 在 FGSCs 和 ESCs 中 拥有 1 或 2 
个 APA 位 点 的 基因 占 最 大 比例 ,总 和 占 总 体 45.2% ,其 
次 是 3 个 和 4 个 ,总 和 为 23.4%, 拥 有 多 至 12 个 位 点 的 基 
因 也 占 10% ,整体 上 看 ,两 者 在 APA 丰 度 上 的 差异 并 不 
显著 。 
2.3 3'UTR 变化 基因 的 功能 分 析 

如 上 文 所 述 , 共 有 795 个 基因 (图 3) 在 FGSCs FŒ 
现 出 相 比 于 ESCs B 3’ UTR 变 短 的 趋势 ,应 用 DAVID 


(https://david.ncifcrf.gov) 的 Gene Ontology 和 KEGG 
对 基因 进行 功能 和 信和 号 通路 分 析 , 这 些 基因 多 富 集 在 和 蛋 
白质 转运 ,定位 , 胞 内 运输 等 生物 学 过 程 ,伴随 着 一 些 豚 
胎 发 育 后 期 出 现 的 与 中 胚层 发 育 和 原 肠 胚 形 成 的 基因 
* 4E RH ;在 细胞 内 的 定位 主要 位 于 胞 内 腔 
(Intracellular lumen) ,细胞 如 内 腔 , 核 内 腔 和 核 质 中 ,以 
及 非 膜 结 合 细 胞 器 ;从 分 子 功能 上 主要 作用 在 与 核 车 酸 
结合 ,主要 是 味 吟 中 的 腺 味 吟 核 背 酸 以 及 一 部 分 与 核糖 
核 昔 酸 和 ATP 的 结合 。 在 通路 分 析 方面 ,这 些 基 因 集中 
于 粘着 链接 ,细胞 周期 以 及 与 卵细胞 减 数 分 裂 相 关 的 信 
号 通路 。GO 分 析 条 目 如 图 4A 所 示 。 
2.4 FGSCs 细胞 内 与 生殖 细胞 发 育 相 关 的 基因 发 生 
3'UTR 缩短 

3 UTR 变 短 的 基因 ,因为 其 较 短 的 3"UTR 增 加 了 
mRNA 本 身 的 稳定 性 ,减少 了 能 与 转录 后 调控 的 反 式 作 
用 因子 (如 RNA 结 合 恒 白 ) 或 者 microRNA 相互 结合 的 
顺 式 作用 元 件 , 更 有 可 能 逃逸 细胞 中 microRNA 介 导 的 
基因 沉默 ,导致 更 高 的 表达 水 平 ” ,对 于 在 FGSCs 中 的 
3' UTR 变 短 伴随 着 表达 量 上 升 的 262 个 基因 ,如 图 3B 
中 所 示 ,在 蛋白 质 转运 和 定位 等 生物 学 过 程 上 有 着 富 
,深入 分 析 发 现在 FGSCs 中 ,与 生殖 细胞 发 育 密切 相 
关 的 基因 fit! 发 生 显著 的 3"UTR 缩 短 。 

Kitl 基 因 编 码 的 KITL 是 一 个 干细胞 因子 (也 被 称 
H SCF) ,作为 配 体 与 跨 膜 酷 氨 酸 激酶 受 体 KIT 结 合 , 两 
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者 形成 kiVkid 信号 系统 共同 调控 着 出 生前 后 卵 梨 中 卵 细 — 3 讨论 
胞 和 卵泡 的 发 育 , 也 有 很 多 报道 在 精子 发 生 中 /8 会 促进 生殖 细胞 的 发 育 是 一 个 复杂 而 高 度 有 序 的 过 程 , 受 


A 型 精 原 细胞 的 DNA 合 成 从 而 促进 细胞 增殖 “2”。 它 在 
FGSCs 中 存在 着 明显 的 3”UTR 缩短 现象 ,在 距离 TITS 
位 点 更 近 的 地 方 存在 APA 位 点 ,如 图 4B 所 示 , 对 近 端 
APA 位 点 的 取 用 伴随 着 该 基因 表达 量 的 升 高 ,表明 jizl 
基因 作为 一 个 参与 生殖 细胞 发 育 的 基因 在 FGSCs 中 有 
特异 的 表达 模式 。 


基因 调控 网 络 精确 和 稳定 的 调控 ” o MERET AUG 
作为 一 种 单 能 性 的 成 体 干细胞 ,处 于 生殖 发 育 进 程 中 一 
个 很 特殊 的 环节 :它们 处 于 原始 生殖 细胞 (PGC) 发 育 下 
游 ,负责 分 化 成 成 的 生殖 细胞 。 近年 来 越 来 越 多 的 研 
究 显 示 APA 在 发 育 过 程 中 发 挥 重 要 的 作用 ,比如 多 能 性 
的 胚胎 干细胞 更 倾向 于 采用 更 短 的 3"UTR ,使 细胞 能 维 
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Fig.4 Gene function analysis of differential 3'UTR genes. A: GO and KEGG pathway analysis using DAVID; B: 
3'UTR shortening in kitl gene shown in IGV genome browse. 


持 自身 的 干细胞 特性 ,保持 自我 更 新 和 较 快 的 增殖 *"。 
而 对 于 生殖 干细胞 的 转录 后 调控 的 研究 仍 处 于 空白 ,在 
此 我 们 利用 基于 高 通 量 测序 的 组 学 研究 方法 获得 了 小 
鼠 FGSCs 在 全 基因 组 范围 内 的 APA 谱 ,使 我 们 能 对 生 
殖 干 细胞 中 的 APA 取 用 模式 有 了 更 全 面 的 认识 ,我 们 
的 工作 填补 了 生殖 细胞 领域 转录 后 调控 机 制 研 究 的 空 
白 ,为 今后 对 某 些 特定 基因 APA 的 深入 研究 提供 了 重 
要 参考 。 

通过 与 小 鼠 ESCs 中 APA 比较 分 析 , 发 现 FGSCs 有 
着 与 ESCs 显著 区 别 的 APA 谱 ,具有 生殖 干细胞 的 特异 
性 :在 我 们 的 数据 中 发 现 有 1 148 个 基因 相对 于 ESCs 发 
Æ T 3 UTR 的 变化 ,其 中 34% 的 基因 存在 3”"UTR 变 长 
现象 ,66% 的 基因 3’UTR 缩短 ,如 图 3A 所 示 。 对 于 3 
UTR 缩短 的 基因 ,在 基因 功能 上 多 富 集 在 与 蛋白 质 转 
运 和 定位 等 生物 学 功能 上 ,这 与 APA 所 介 导 的 与 蛋白 
质 细 胞 定位 相关 功能 相 一 致 ,在 一 些 豚 胎 发 育 后 期 出 现 
的 与 中 胚层 发 育 相 关 基 因 条 目 也 有 高 集 ( 图 4A) ,表明 
这 些 在 豚 胎 发 育 阶段 发 挥 作 用 的 基因 可 能 在 生殖 干 细 


胞 里 也 发 挥 一 定 的 作用 。 

AF f] 3 UTR 能 导致 更 高 的 转录 活性 ,其 转录 激 
活 是 通过 顺 式 作用 元 件 的 缺失 ,逃逸 microRNA 介 导 的 
负 调 控 来 实现 中 。 针 对 上 肉 性 生殖 干细胞 中 3 UTR 缩短 
的 基因 我 们 展开 进一步 分 析 , 发 现 其 中 包含 一 些 与 生殖 
发 育 密切 相关 的 基因 ,比如 fit 基因 (图 4B), 该 基因 编码 
蛋白 作为 开 型 跨 膜 酷 氨 酸 激酶 受 体 KIT 和 蛋白 的 配 体 "”， 
A BE Sp E TEE SERT IV dE 
FGSCs rf ,该 基因 产生 了 APA 介 导 的 37UTR 缩短 转录 
子 ,而 且 表 达 量 较 高 ,推测 缩短 的 3"UTR 导致 该 基因 的 
转录 激活 而 参与 生殖 发 育 相 关 的 生物 学 活动 。 

尽管 3`UTR 缩短 导致 转录 激活 ,我 们 在 FGSCs 中 
EERI UTR 变 短 的 基因 并 不 总 是 伴随 着 高 表达 量 : 
7954^- 3 UTR 缩短 的 基因 中 ,其 表达 量 较 ESCs 上 升 的 
只 有 262 个 ,剩余 的 533 个 基因 没有 明显 上 调 甚 至 表现 
出 表达 下 调 的 趋势 ,如 图 3B 所 示 , 这 提示 3’UTR 的 长 
短 与 表达 量 不 仅仅 是 简单 的 负 相 关 。 除 了 我 们 的 结 
外 ,之 前 的 一 些 报 道 也 发 现 类 似 的 现象 “1, 我 们 推测 
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对 于 那些 3’UTR 缩短 但 是 表达 量 没 有 升 高 的 基因 来 
说 ,可 能 是 细胞 内 的 microRNA 没 有 表达 ,从 而 不 能 有 效 
的 通过 与 3"UTR 和 靶 区 域 结合 抑制 靶 基 因 的 表达 所 致 。 

综 上 所 述 ,本文 首次 详细 摘 绘 了 小 鼠 FGSCs 在 全 
基因 组 范围 内 的 APA 谱 ,通过 和 小 鼠 胚胎 干细胞 进行 
系统 性 的 比较 ,为 从 转录 后 调控 水 平 更 深入 地 了 解 生殖 
干细胞 的 基因 调控 机 制 提供 了 基础 ,从 而 进一步 丰富 了 
生殖 干细胞 生物 学 研究 内 容 。 
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